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We study the transport properties of a graphene ferromagnet-insulator-superconductor (FIS) junc- 
tion within the Blonder-Tinkham-Klapwijk formalism by solving spin-polarized Dirac-Bogoliubov- 
de-Gennes equation. In particular, we calculate the spin-polarization of tunneling current at the 
I-S interface and invesigate how the exchange splitting of the Dirac fermion bands influences the 
characteristic conductance oscillation of the graphene junctions. We find that the retro and spec- 
ular Andreev reflections in the graphene FIS junction are drastically modified in the presence of 
exchange interaction and that the spin-polarization (Pt) of tunneling current can be tuned from 
the positive to negative value by bias voltage (V). In the thin-barrier limit, the conductance G of a 
graphene FIS junction oscillates as a function of barrier strength \. Both the amplitude and phase 
of the conductance oscillation varies with the exchange energy E ex . For E ex < Ep (Fermi energy), 
the amplitude of oscillation decreases with E ex . For E% x > E ex > Ep-, the amplitude of oscillation 
increases with E ex , where E%. = 2Ep + Uo (Uo is the applied electrostatic potential on the super- 
conducting segment of the junction). For E ex > E% x , the amplitude of oscillation decreases with 
E ex again. Interestingly, a universal phase difference of n/2 in \ exists between the G — \ curves for 
E ex > Ep and E ex < Ep. Finally, we find that the transitions between retro and specular Andreev 
reflections occur at eV = \Ep — E ex \ and eV = E ex + Ep, and hence the singular behavior of the 
conductance near these bias voltages results from the difference in transport properties between 
specular and retro Andreev reflections. 

PACS numbers: 74.45.+C, 72.25.Mk, 73.40.Gk, 73.63.-b 



I. INTRODUCTION 

Graphene is a single atomic layer of graphite with 
a two-dimensional (2D) honeycomb lattice structure [l|. 
The valence and conduction bands of graphene touch 
each other at six corner points of the 2D hexagonal Bril- 
louin zone. These six corner points are known as Dirac 
points and are divided into two incquivalcnt groups, de- 
noted by K and K', respectively. In 1947, Wallace^ 
first predicted that graphene possesses a linear Dirac- 
like energy dispersion near each Dirac point. The lin- 
ear approximation for energy dispersion is valid even up 
to 1 eV above and below the Dirac points. However, 
the fabrication of graphene was realized only a few years 
agoQ. The linear dispersion relation of graphene has 
been recently observed in angle-resolved photocmission 
spectroscopy (ARPES) experiment Q ■ Due to the pe- 
culiar electronic structure, charge carriers in graphene 
behave like massless relativistic particles (called Dirac- 
like fermions) at low energes. Recent theoretical works 
showed that Dirac-likc fermions in graphene may result 
in several novel phenomena such as anomalous quantum 
Hall effect @ and Klein tunnclinglg, Ml. The anomalous 
quantum Hall effect [1, Q and Klein tunneling [ToL [Til. [l2j 
have been observed recently. In addition, the Fermi level 
of graphene can be tuned by a gate voltage through the 
electric field effect Q. It is also possible to create fer- 
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romagnctism and superconductivity in graphene via the 
proximity effect. 

Superconductivity is one of the interesting proper- 
ties of graphene. In nature, graphene is a zero- 
gap semiconductor and not a superconductor. How- 
ever, superconductivity of graphene may be induced 
by placing a superconducting electrode on the top of 
graphene via the proximity effect [H, EH EE Htl Il7l 
or coating a sheet of graphene with metallic atoms |18j. 
There have been many studies on superconductivity of 
graphene. Beenakkerpjj showed that a normal metal- 
superconductor (NS) graphene iunction[2fj| can exhibit 
unique specular Andreev reflection which is qualitatively 
different from usual Andreev reflection in a conventional 
NS junction [2l]]. Josephson current of a superconductor- 
normal metal-superconductor (SNS) graphene junction 
has also been calculated [13, and observed [24]. The 
observation of Josephson current confirms that supercon- 
ductivity of graphene can be induced via the proximity 
effect. It was shown that the tunneling conductance of 
a normal mctal-insulator-superconductor (NIS) graphene 
junction oscillates with barrier strengthfcalled the con- 
ductance oscillation effects) [IE HEIIS 28]. This behavior 
is analogous to the Klein tunneling. Notice that the insu- 
lator here is not an insulator with an energy gap between 
the valence and conduction bands. The notation " I" in 
a graphene NIS junction means just a normal segment 
of graphene with a large electrostatic potential. Such 
an insulator could be created by using a gate voltage or 
chemical doping [||. 

Recent theoretical studies showed that intrinsic for- 



romagnctism of grapheme may exist [29l. 130 . l3l| . but 
it has not been observed so far. Room-temperature 
ferromagnetism was observed [HJ, but it was attribued 
to the presence of the defects on graphene. However, 
proximity-induced ferromagnetism in graphene was re- 
cently realized experimentally [33j| . Furthermore, a theo- 
rectical work showed that the spin injection efficiency for 
C0/AI2O3 /graphene could amount to 18% and could also 
be increased up to 31% by applying a current bias(34|. 
Therefore, it is feasible that a segment of graphene could 
be ferromagnetic. Recently, quantum coherence trans- 
port in a graphene ferromagnet-superconductor (FS) 
junction was investigated theoretically [35l 13(1 \3l\ . These 
studies suggested that the conductance of graphene could 
increase with the exchange energy E ex when the E ex is 
larger than the Fermi energy Ep . This finding is dramati- 
cally different from conventional FS j unctions [H|,[3!| and 
strongly contradicts our intuition. As stated above, the 
interplay between superconductivity and ferromagnetism 
of graphene is interesting. The tunneling behavior of 
Dirac-like quasipariticles, the analogous Klein tunneling, 
is also remarkable. 

In this paper, therefore, we study the transport proper- 
ties of a graphene fcrromagnct-insulator-supcrconductor 
(FIS) junction by using Blonder-Tinkham-Klapwijk 
(BTK) formalism [2l|. The previous studies on the trans- 
port properties of graphene F/S juctions focus on the ef- 
fect of the exchange energy on zero-bias conductance [35| 
and conductance spectra 36]. In contrast, here we calcu- 
late the spin-polarization of tunneling current at the I-S 
interface and also invesigate how the exchange splitting of 
the Dirac fcrmion bands influences the characteristic con- 
ductance oscillation of the graphene junctions. In partic- 
ular, we find that the spin polarization of a graphene fer- 
romagnetic segment is not always positively correlated 
with the exchange energy and that the tunneling spin 
polarization could be negative. In addition, we calculate 
the conductance in the thin-barrier limit, where thick- 
ness d — > and the barrier potential Vq — > 00. In this 
limit, the conductance (G) of a FIS junction oscillates as 
a function of the barrier strength x an d the amplitude of 
oscillation of conductance varies with E ex and is minimal 
at E ex = Ep. The G — x curve for E ex > Ep is phase 
shifted by tt/2 in barrier strength %, compared with the 
curve for E eT < Ep. 



II. THEORETICAL METHOD AND 
ANALYTICAL CALCULATIONS 

A. Theoretical model 

Here we consider a graphene FIS junction, as shown 
in Fig. 1(a), with region F in x < —d, region I in 
—d<x < and region S in x > 0. To tackle a problem 
concerning superconductivity and relativity, we should 
consider the Dirac-Bogoliubov-de-Genncs (DBdG) equa- 
tion. The DBdG equation for the graphene nonmagnetic 
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FIG. 1: (Color online) (a) A schematic plot of a graphene FIS 
junction. Here, F, I, and S denote the ferromagnetic, insulat- 
ing and superconducting regions, respectively. The F and S 
regions could be manufactured by contacting graphene with 
a ferromagnetic electrode and a superconducting electrode, 
respectively. The I region could be created by an external 
gate voltage Vb- An additional gate voltage Uo may be ap- 
plied on the S region, (b) A schematic plot of the electrostatic 
potential U(r) in the FIS junction of (a). 



metal-superconductor (NS) interface was derivied in Ref. 
fl9l ] by using a time-reversal operator T = (gv<%>7v)C|40l]. 
Here C is the complex-conjugation operator, aj and Tj 
are the 2x2 Pauli matrices acting on sublattices (i.e. 
pesdospins) and valleys, respectively. The final form of 
the DBdG equation in [lj| was decomposed into two de- 
coupled four-dimensional sets 



H± — Ep\ 

A*i 



Al 

(H ± - Epl) 



(1) 

where A denotes the pair potential in the superconduct- 
ing segment and 1 is the 2x2 indentity matrix. H± is 
the single particle Hamiltonian with indices ± labeling 
the two valleys of the band structure at K and K', re- 
spectively. The single particle Hamiltonian H± is given 

byQ 



H± = -ihv F (a x d x ± dydy) + U. 



(2) 



where vf represents the Fermi velocity with v F m 
10 6 m/s[l|. In Eq. (1), u and v denote the wave func- 
tions of the electron- like and hole-like excitations, respec- 
tively. It is obvious that both u and v are composed of 
two components which represent sublattices A and B in 
grephene, respectively. Note that the solutions of Eq. 
(1) must satisfy the condition: (y~ v + ) T = T{u + u~) T , 



3 



where 



") T and 



are, respectively, the so- 



lutions of Eq. (1) with H + 



and In other words, 

v~ = a z Cu~ and v + = a z Cu + . As can be seen from Eq. 
(1), the electron-like part of wavefunction of K valley u + 
is coupled with v~, which is related to the electron-like 
part of wavefunction of K' valley u~ , by the energy gap 
A. A similar relation exsists between u~ and v + . This 
indicates that the two valleys K and K' are coupled to- 
gether by A for a graphene superconducting segment. 

To take ferromagnetism into consideration, we regard 
the ferroma gne tic segment as a free-electron-like Stoner 
ferromagnet [4l| . Its spin-up band is shifted down by the 
exchange energy h, while its spin-down band is shifted 
up by h. The intrinsic spin-orbit coupling in graphene is 
negligible [53]. Following Ref. [l9| and introducing real 
spin degrees of freedom, we find that the DBdG equation 
in the presence of an exchange interaction is given by 



H± - E F 1 - p a hl Al 

A*i -(h±- E F i + Pa hi) 



(3) 



Subscript cr(a) labels the real spin. If a denotes up 
(down) spin, a denotes down (up) spin, pj is 1 and pj 
is -1. Note that the indices |, J, refer to the two spin 
subbands, respectively. In addition to the coupling be- 
tween the two valleys, Eq. (3) also implies that an An- 
dreev reflected hole is formed due to the removal of an 
electron whose spin direction is opposite to that of the 
incident electron. The hole is an antiparticle of the re- 
moved electron. Therefore, the hole has the same spin 
direction as the incident electron. Let us call Eq. (3) 
the spin-polarized DBdG equation and the details of its 
derivation are given in Appendix A. 

We assume that ferromagnetism and superconductiv- 
ity are uniformly induced in region F and region S, re- 
spectively. Therefore, h, A, and U are modeled as step 
functions, 



h{v) = E ex Q{-x-d), 



A(r) = A e^e(a 



U(r) = -U Q{x) + V Q(x + d)Q(-x), 



(4) 



where Q is the Heaviside step function and <f> is the 
phase of the pair potential of the superconducting seg- 
ment. The potential U is plotted in Fig. 1(b). How- 
ever, the DBdG e quat ion is derived within the mean- 
field approximation [43j. The mean-field approximation 
is valid only when the Fermi wave length A' F in region S is 
much smaller than the superconducting coherence length 
f, where £ cx 1/AoH and X' F cx 1/(E F + U ). There- 
fore, our calculations should be done under the condition: 
Aq <C (E F + Uq), When the order of magnitude of E F 
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FIG. 2: (Color online) The excitation spectra in the ferro- 
magnetic graphene segment in (a) for E ex < Ef and in (b) 
for E ex > Ef, calculated using Eq. (5). Red lines indi- 
cate electron excitations, while blue lines indicate hole exci- 
tations. Solid and dotted lines denote the conduction and 
valence bands, respectively. 



O(Ef) < O(Ao), wc must adopt the regime Uo 3> Ao. 
We can always satisfy the requirement of the mean field 
approximation by adjusting Uq via the gate voltage. 

Because of the valley degeneracy, we only need to con- 
sider the sets with H + in Eq. (3). We solve Eq. (3) with 
H+ and obtain the dispersion relations of excitations in 
region F 

e eT = hv F \k.\ - E ex - E F , 



'-hi 



±hv F \k\- E ex + E F , 



= ±hv F \k\ +E ex 



E F , 



±Hv F \k\+ E ex + E F , 



(5) 



where subscript e (h) indicates electron (hole) excita- 
tions. + and — represent the dispersions of the bands 
plotted with dotted lines and solid lines, respectively, in 
Fig. 2. The crossing points between the conduction and 
valence bands are at e = \E F — E ex \ and e = E F + E ex . 

We also obtain the dispersion relation of quasiparticles 
in region S 



e s = J(E F + U ± hv F \k\) 2 + Ag. 



(6) 



In addition, the wave functions in region F, region I and 
region S are also obtained. 



B. Calculation of conductance 

We can calculate the Andreev reflection (a a ) and nor- 
mal reflection (b a ) amplitudes by applying the appropri- 
ate boundary conditions to match the wave functions. 



The details of the calculation are described in Appendix 
B. The tunneling conductance is calculated by consid- 
ering the contributions of the incident spin-up and spin- 
down electrons. The zero-temperature tunneling conduc- 
tance of the FIS junction is given by 

G(eV) = ^G a {eV) I ' da a cos{a a )x 



(a) modified specular Andreev reflection (b) modified retro-Andreev reflection 



cos(a CT ) 



2e 2 

G a (eV) = —N„{eV), N a {e) 



\b a (eV,a„)\ z , 



\e + E F + p a E ex \W 

HvpTT 

(7) 

where G a is the normal-state conductance of the in- 
cident spin-er electron. N a is the spin-cr density of 
state (DOS) in a graphene sheet of width W. Finally, 
we calculate the normalized conductance G/Gq, where 
Go = J2a Ga = Gf + Gi- However, a critical angle of 
incidence a c a exist in Eq. (7) because the wave functions 
would decay as the incident angle a a > a.%. Adopting 
an approach used in Ref. (HI, we do not consider the 
evanescent wave functions and hence set the integration 
limit to the maximum angle of incidence in Eq. (7). The 
critical angle can be written as 

o£(e) = min« 1 (e),< 2 (e),< 3 (e)), 



0) = sin 1 { 



V3 



(e) = sin- 1 ( 





- P*E ei 


c - E F \ 
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-E F + 


P*jE ex | 


\e 


- \E F 


-Veil 




-E F + 


PcrE ex 




(E F + 


U ) 
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-E F + 


PaE ex 



(8) 



When a a > a c all the wave function of the Andreev re- 



flected hole decays. For a a > 



V2! 



the quasiparticlc 



in the barrier region cannot propagate. For a a > a£ 3 , 
the wave functions of the transmitted quasiparticles are 
evanescent. If sin(a^ 1 ) is larger than 1 by using Eq. (8), 
the Andreev reflected hole do not decay at any angle of 
incidence. Similarly, for sin(a£. 2 ) > 1 [ sm ( a cr3) > l]i the 
quasiparticle can always propagate in region I [region S] . 
Since we do not take any decaying state into account, the 
critical angle is the minimum among a% 1 , OL c a2 , and a^ 3 . 



C. Modified specular and retro-Andreev reflection 

We note that in the Andreev reflection process, the 
emerging hole is usually reflected back along the path of 
the incoming electron [45j. Nevertheless, this may happen 
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|kha|. If |k ecr | is not equal Otherwise, v e 



FIG. 3: (Color online) Schematic plots: (a) the modified spec- 
ular Andreev reflection, (b) the modified retro-Andreev reflec- 
tion, and fc) the normal refletion. 



to |kho-|, the Andreev reflected hole does not retrace the 
path of the incident electron anymore Hi . However, for 
a graphene superconducting junction, the hole dose not 
neccessarily go along the path of the incident electron 
even if |k ecr | = |kha|- Two different types of Andreev re- 
flection may happen. One is the Andreev retroreflection 
and the other is the specular Andreev reflection In 
the heavily doped graphene (with E F 3> A ), the An- 
dreev retroreflection takes place and the outgoing hole 
retraces the path of the incident electron. In the weakly 
doped graphene (with E F -c A ), the specular Andreev 
reflection takes place and the outgoing hole travels like 
the normal reflected electron. The specular Andreev re- 
flection is a unique phenomenon in a graphene supercon- 
ducting junction. 

For a two or three dimensional system, the boundary 
conditions would require that there is no scattering in 
the transverse directions (i.e., the directions parallel to 
the interface). That is to say, the transverse wave vec- 
tor is conserved in all scattering processes. For the sys- 
tem investigated here, the transverse direction means the 
y-dircction. The wave vectors of the incident electron 
and scattered particles in the y-dircction are all equal 
to q, as shown in Appendix B. From this, we deduce 
that |k ecr sm(a&) \ = \khs sin(a s ) |. The incident angle 
a a equates to the Andreev reflected angle ot s only when 
|ketr| = |khff|- Using v = Vke/ft and Eq. (5), where v 
denotes the group velocity, we can determine which di- 
rection an excitation goes. If (e ecr + p a E ex + E F ) > 0, 
e ea = hv F \k\ — p a E ex — E F and the velocity of the in- 
cident electron v e is parallel to k. Conversely, v e is an- 
tiparallel to k. Similarly, if {s] lS + p a E ex — Ep) > 0, 
£fis = hv F \k\ — p a E ex + Ep and the velocity of the re- 
flected hole Vh is parallel to k. Conversely, Vh is an- 
tiparallel to k. Therefore, if both (e ecr + p a E ex + Ep) 
and (sua + PaEex ^ Ep) have the same sign, v ey = Vh y - 
—Vhy The relation between a a and a s 



5 



is given by 



III. SPIN POLARIZATION 



. _! (s + p a E, 
sin ( 



Ep) 



(e + p a E ex - E F ) 



(9) 



Here, the incident angle a a is always positive while a a 
may be either positive or negative. If a s is positive, 
v ey = Vhy and the modified specular Andreev reflection 
happens. If c% is negative, v ey = —V} ly and the modified 
retro- Andreev reflection occurs. Here, " modified" means 
that a s ^ a a . 

Kashiwaya et al. reported that under the effect of 
the exchange energy the Andreev reflection angle is not 
equal to the incident angle any more(46j. Interestingly, 
for a doped nonmagnetic graphene superconducting junc- 
tion with e comparable to Ep and E ex = 0, a s is not 
equal to a a either. The Fermi energy of a conventinal 
metal, such as Cu or Al, is in the order of I eV(47| 
and hence is several thousand times as big as the bias 
voltage e. For conventional nonmagnetic metals, there- 
fore, |k etT | »2 |kha| ~ kp, where the Fermi wave vector 
kp = \/2mEp/h, and the incident angle a a is always 
equal to the Andreev reflected angle a s . It is interest- 
ing that for a graphene superconducting junction, when 
e~Ep — p a E ex , the path of the Andreev reflected hole 
would severely deviate from the normal line at the small 
incident angle. By changing the Fermi energy via a gate 
voltage, we can reach the condition of e~Ep — p a E ex eas- 
ily. Fig. 3 illustrates the modified retro- Andreev reflec- 
tion, specular Andreev reflection and normal reflection. 
From Fig. 3, it can be seen that the normal reflection in 
a graphene superconducting junction is not qualitatively 
different from a conventional superconducting junction. 
However, the retro- and specular reflectivity of Andreev 
reflection can be drastically broken under the effect of 
the exchange energy and the bias for the graphene su- 
perconducting junction. 

In addition, when either (e + p a E ex — Ep) or (e + 
p a E ex + Ef) changes sign, there is a trasition between 
the retro- Andrcevc reflection and specular Andreev re- 
flection. This happens at a Dirac point. Take Fig. 
2(b) for example, the Dirac points of the bands of the 
spin-down electron and hole excitations are located at 
e = (E ex ~Ep) and e = (E ex +Ep), respectively. An inci- 
dent spin-up electron proceeds with the specular Andreev 
reflection when e < (E ex — Ep)- As the bias voltage e is 
increased up to the Dirac point of the band of the spin- 
down electron (E ex — Ep), the electron would undergo 
the retro- Andreev reflection. Keeping increasing the bias 
voltage, the electron would proceed with the specular An- 
dreev reflection again when e > (E ex + Ep). Similarly, a 
trasition between the retro- Andreeve reflection and spec- 
ular Andreev reflection occurs at the Dirac points as a 
spin-up electron incidents. In short, for a graphene FIS 
junction, there are two transition points. This is quite 
different from a graphene NIS junction which has only a 
transition point. 



A. Spin polarization of DOS 

Measurement of spin polarization is an important sub- 
ject in spintronics. There are several ways to mea- 
sure spin polarization, such as spin-polarized photoemis- 
sion spectros cop y pjj . Meservey-Tedrow spin tunnelling 
spectroscopy [43 and Andreev reflection spectrum [50|. 
The BTK formula provide a quite simple method to 
measure the spin polarization of a metal by Andreev 
reflection spectrum. How to define the spin polariza- 
tion is also a valuable question and has been discussed 
before by Mazin[5lJ. The most popular definition of 
spin polarization in a ferromagnet is P = (Ni(Ef) — 
N l (EF))/{N r (E F ) + N^Ep)), where N^Ep) [N L (E F )} 
denotes spin-up [down] DOS at the Fermi level. For a 
conventional metal, the Fermi level does not move much 
when a bias is applied, and the DOS remains almost con- 
stant. Now consider the effect of the bias on the DOS for 
graphene. iVj and N± are now given in Eq. (7). There- 
fore, the spin polarization of DOS can be written as, 



P 



{e + E F + E e 



\e + E F - E e 



(e + E F + E e 



\e + E F - E e 



(10) 



The spin polarization in a graphene ferromagnetic seg- 
ment varies with the applied voltage. When e < E ex — 
E F , P = (e + E F )/E ex . Conversely, P = E ex /(e + E F ). 
In the low bias regime (E cx ,Ef ^S> e), we find that 
P = E ex /E F HE ex < E F and P = E F /E ex if E ex > Ep. 
The spin polarization becomes independent of s. Fur- 
thermore, it is counterintuitive that the spin polarization 
is not always positively correlated with E ex . It is clear 
that we could increase the spin polarization in graphene 
by tuning the Fermi energy via a gate voltage. The Fermi 
energy in graphene is usually in the order of 0-100 meV. 
For the lightly doped graphene, therefore, we only need 
to apply a gate voltage of a few milli-clectron voltage 
to change the spin polarization in graphene significantly. 
Furthermore, for Ep — E ex 3> £, graphene is a half metal 
with P = 1. The exchage energy is fixed by the proxim- 
ity effect. Hopefully, this interesting prediction of half- 
metallic graphene will stimulate future experiments on 
graphene by tuning the gate voltage to the condition of 
Ef = E ex . 

Nonetheless, it is worth mentioning that when Ep is in 
the vicinity of the neutrality point (i.e. the Dirac point), 
the carrier density goes to zero. The electron- hole pud- 
dles can appear at low carrier densities and in the pres- 
ence of disorder. [H, HH In this case, local variations of 
the Fermi energy should be taken into account. It is 
well known that a perfect graphene sheet is very difficult 
to obtain experimentally 53]. Therefore, in experiments, 
corrugations are one source of disorder. Fascinatingly, 
it has been very recently demonstrated that an ultra- 
flat graphene with height variation less than 25 pm can 
be obtained [H3], suggesting that disorder induced behav- 
iors such as electron-hole puddle formation, of graphene 
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FIG. 4: (Color online) (a)The tunneling spin polarization for 
several values of the exchange energy E ex with Ef = lOOAo, 
Uo = and \ = t/2. The insert: the Px-eV curve with the 
same parameters as (a) and eV in the range of 0-500Ao. For 
E ex < Ef, the spin polarization could be negative. (b)The 
effect of barrier on the curve of the tunneling spin polarization 
versus bias voltage for EF,E ex ^> Aq with Ef = lOOAo, 
E ex = 50A , and U = 0. 



could be avoided. In a ferromagnetic graphene, the car- 
rier density of spin-down electrons would approach zero 
in the vicinity of E ex at Ep = E ex . However, unlike in 
nonmagnetic graphenes, the Dirac points of different spin 
species in a ferromagnetic graphene are not located at the 
same position. Therefore, the total carrier density is not 
necessarily close to zero as Ep approaches E ex , and our 
calculations would be still valid in the sufficiently clean 
(little dirty) regime. 



B. Tunneling spin polarization 

The spin polarization of tunneling current for 
graphene FIS junction is defined as 



Pi 



(*T ~ h) 



(ii) 



where I a is the spin-cr current which injects into the su- 
perconductor at the I/S interface. The detailed calcu- 
lation of Pt is displayed in Appendix C. It should be 
emphasized that the magnitude of the tunneling spin po- 
larization Pt here represents the spin injection efficiency 



of a graphene FIS junction at the I/S interface. Due to 
spin diffusion [55t . the spin current in the junction docs 
not remain constant, and indeed, deep inside the super- 
conductor, Pt of the tunneling current would go to zero. 

We show the spin polarization of tunneling current as 
a function of bias voltage in Fig. 4 (a). Here, we pay 
more attention on the bias in the range of — 30Ao be- 
cause the Andreev reflection spectrum is usually mea- 
sured at eV w — 10 meVpSfl [5tjj . The enegy gap A 



is around 0.5 meV|26[. Nonetheless, Pr-eV curves with 
eV in the range of — 500 meV are also plotted in the 
insert in Fig. 4(a). From Fig. 4(a), we find the tun- 
neling spin polarization Pt is always zero at subgap bias 
for any E ex . This is because only supercurrent can flow 
through the interface for E < Ao and singlet supercur- 
rent carries no spin polarization. For E ex < Ep, exclud- 
ing E ex fa Ep, Pp slightly increases and then decreases 
with increasing eV up to around eV = Ep — E ex . Tak- 
ing the curve with E ex = 50Ao as example, at around 
eV — Ep — E ex , the tunneling spin polarization reaches 
its minimum. Futhermore, after Pt reaches the minimum 
at around eV = Ep — E ex , Pt is raised by eV again. Un- 
till around eV = Ep + E ex , a local maximum appears, 
and the Pt then decreases monotonically as eV increases. 
Finally, Pt approaches + at high bias. Interestingly, Pt 
would become negative when the bias is increased up to 
some critical value M (Fig. 4a), which is the transi- 
tion point between the positive and negative tunneling 
spin polarization. However, for E ex > Ep, excluding 
E ex Ep, Pp increases up to its maximum and then 
falls. The maximum happens at around eV = E ex — Ep. 
After Pp falls to some degree, it would increase with 
eV again untill eV is about eV = E ex + Ep. A local 
maximum occurs at about eV = E ex + Ep. After Pp 
reaches the local maximum at around eV = E ex + Ep, 
it decreases monotonically with increasing eV. It is clear 
from Fig. 4(a) that unlike the case with E ex < Ep, the 
tunneling spin polarization for E ex > Ep is always posi- 
tive. 

In short, the magnitude of tunneling spin polarization 
\Pt\ reaches its maximum at around eV — \Ep — E ex \ 
and a local maximum at around eV = Ep + E ex for 
any E ex which is not close to Ep and the graphene FIS 
junction can exhibit the unique negative-value tunneling 
spin polarization if E ex < Ep. The negative tunnel- 
ing spin polarization implies that it is harder for major- 
ity carriers (spin-up electrons) to tunnel through the in- 
terface than for minority carriers (spin-down electrons). 
eV = \Ep — E ex \ and eV = Ep + E ex are the Dirac 
points of the spin-down electron and hole bands, as shown 
in Fig. 2. The Andreev reflection is suppressed at the 
Dirac points of the hole bands. No spin-polarized cur- 
rent can flow through the interface via the Andreev re- 
flection process. Therefore, near the Dirac points of the 
hole band, spin-polarized current is easier to be carried 
over and hence the tunneling spin polarization can be 
large when eV is located at the Dirac points of the hole 
bands. Near the Dirac points of the spin-down electrons, 
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FIG. 5: (Color online) (a) The conductance spectra for sev- 
eral values of the exchange energy E ex with Ef — lOOAo, 
Uo — and \ — n/2. (b) The conductance spectra for sev- 
eral values of barrier strength \ with Ef = 100 Ao, E ex — 
50Ao, and Uo = 0. The curves for x = and x — T 
are identical. Note that Go is not always constant. For 
Ef, E ex ^>eV, Go is independent of eV. When E ex < Ef, 
Go = 4(£ F /A )(e 2 W//i.Co), where £ = irHv F /A . Con- 
versely, Go = 4(E ex /Ao)(e 2 W/hto). 



the spin-down DOS is very small and therefore spin-up 
current predominates. In this case, the spin-down cur- 
rent is small and hence the tunneling spin polarization 
would be also large. We find that the bias voltage can 
be used to tunc the tunneling spin polarization. Taking 
the case with E ex = 120 An. as an example, Pp can vary 
from 0% to 62% (i.e., - 0.62) in the range of - 30 
meV. The effect of barrier on the magnitude of tunneling 
spin polarization is shown in Fig. 4(b). The magnitude 
of the spin polarization change under the effect of bar- 
rier strength can be, e.g., 9% (0.09) at eV = 30A for 
the case with Ep = 100A , E ex = 50A , and Uo = 0. 
Furthermore, we notice that x could move the position 
of the transition point M (Fig. 4b). 



IV. NUMERICAL RESULTS 
A. Conductance in the large Fermi energy limit 

In the limit of a thin barrier, where d — * and Vq — ► oo, 
we use Eq. (7) and Eq. (A9) to calculate the conduc- 
tance. First, we consider the case of Ep ^> Ao. Because 
Vb 



oo, a c a2 must be larger than a c al and o^ 3 , where 



V2 

and 



V3 



arc defined in Sec. II. Here, we are 
interested in the conductance for a small bias voltage 
(0 < e < 2Ao). Using Eq. (8), we find that is always 
equal to 90° in the low bias regime. Therefore, we only 
need to discuss af. If E ex < 2Ep 

We 

Uo). 
The 

5-7. 



r3 



\-U , 
and 



cv5 3 and 



Q'l = o^j. Conversely, > a' 
define the critical exchange energy E^ x = (2Ep 
Here, we adopt Ep = 100A > A and U = 0. 
critical exchange energy is equal to 200 Ao. 

The results with Ep 3> Ao are displayed in Figs. 
As a check of the validity of our calculations, we empha- 
size that the results of a graphene NIS junction reported 
in Ref. [28| can be reproduced by setting E ex = 0. It 
can be seen by comparing the curve for E ex = in Fig. 
5(a) with the results of Ref. [28[. We can also reproduce 
the previous results of a graphene FS junction from Ref. 
[HI and Ref. [36| by setting x = 0- This also indicates 
that our formalism is more general. From Fig. 5(a), 
we find that for the large Ep (Ep Ao), the conduc- 
tance of a graphene FIS junction increases with E ex for 
E c ex > E ex > E F . Conversely, for E ex < Ep, the con- 
ductance would decrease with increasing E ex . However, 
for E ex > E^ x , the conductance decreases with increas- 
ing E ex again. In principle, the spin polarization (P) of 
DOS would suppress the Andreev reflection [4l|. How- 
ever, when Ep e (i.e., e — > 0), Eq. (10) indicates that 
P increases with increasing E ex for E ex < Ep and de- 
creases with increasing E ex for E ex > Ep if Ep is fixed. 
Therefore, when E ex < E c ex , the conductance would show 
the behavior as seen in Fig. 5(a). At E ex = Ep, P 
reaches its maximum value of 1.0, and hence the con- 
ductance becomes very small. However, for E ex > E° x , 
although E ex does not lower the probability of Andreev 
reflection, it would obstruct the propagation of quasipar- 
ticles in the superconducting segment. As E ex > E° x , a*j 
is equal to a| 3 . a| 3 decreases with increasing E ex and 
hence the transmission of quasiparticles in region S falls 
as E ex raises. Therefore, E ex reduces the conductance 
even if P is enhanced by E ex . 



B. Conductance vs barrier strength 

From Fig. 5(b), we find the conductance still oscil- 
lates as a function of x with a period tt like a graphene 
NIS junction under the effect of exchange energy. This 
can also be understood from Eq. (A9). Both the nor- 
mal and Andreev reflection amplitudes (a a and b a ) are 
oscillatory functions of x with a period 7t. Therefore, 
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FIG. 6: (Color online) The zero-bias conductance as a func- 
tion of barrier strength x with Ef = 100 Ao, Uo = 0, and 
E ex = 0, 50A , IOOAq, 200A , 300A . 



the conductance G(eV) oscillates as a function of x with 
a period 7r. However, this unique oscillation originates 
from the relativity of quasiparticles in graphene. Here 
the relativity means that quasiparticles in graphene pos- 
sess a linear energy dispersion. Our calculation shows 
that the relativity is not destroyed by exchange split- 
ting. Dirac-like quasiparticles can transmit through a 
high barrier, in contrast to nonrelativistic particles. The 
motion of non-relativistic particles is described by the 
Schrodingcr equation. Adopting the step-barrier model 
used in our calculation (see Fig. 1(b)), we find the wave 
functions of nonrelativistic particles always decay in the 
barrier region for E < Vq. The wave functions of non- 
relativistic particles are traveling waves and resonance 
scattering [53] may happen only for E > Vq. However, 
for Dirac-like quasiparticles, the wave functions are trav- 
eling waves and resonance scattering may happen even if 
E < Vq. As the wave functions are traveling waves, the 
transmission coefficient must be an oscillatory function of 
kbd, where kb is the wave vector of particles in the barrier 
region. We let Vq — > oo and obtain kf,d — > \- Therefore, 
the conductance oscillates with \ even if Vq — > oo. 

Next, let us consider the zero-bias conductance as a 
function of barrier strength \, plotted in Fig. 6. The 
zero-bias conductance oscillates with \ over a period of 7r. 
The oscillatory behavior results from the Klein tunneling 
in a superconducting junction. The Klein tunneling has 
two remarkable characteristics: i) the transimission is 
not always suppressed by a barrier and ii) for normal 
incidence, the barrier is perfectly transparent. Both of 
the two characteristics have been observed [l^, EH- I n 
addition, we find that for E ex < Ep, G(0) reaches its 
maximum at x = n7: where n is an integer (i.e., n — 
0,1,2,...). In contrast, for E ex > Ef, G(0) reaches its 
maximum at \ = ( n + 1/2)tt, i.e., G(0) is phase shifted 
by 7r/2 in x, compared with that for E ex < Ep. Due 
to the phase shift of 7t/2, in the reange of x = 0-7r/2, 
for E ex > Ep, the conductance increases steadily with 



2.5 
2 

1.5 
1 

0.5 


2 

1.5 
1 

0.5 



Ja) 

E F =100A /.* .X 



e=0 



/ 

/ .* 



— x=o 

JC=ji/4 

■ ■ %=n!2 



tii- 



* \ 

\ \ 
♦ \ 
♦ \ 



-50 



50 



~^ 1 

(b) 

E =50A A 

ex 

E c =100A n 

F 

e=0 
U o=° 



T 




— x=o 

JC=ji/4 

■ ■ X=t/2 



50 



/ ex . \ 

/ . E =100A n ♦ \ 

F 



— x=o 

— - x=k/4 
■ ■ X=%ll 






a( ) 



50 



FIG. 7: (Color online) (a) The transmission coefficient as a 
function of incident angle a in a NIS junction with e = 0, 
Ef = lOOAo, and Uo = for \ — 0, 7r/4, tt/2, respectively, 
(b) The spin-up transmission coefficientin Tf in a FIS junction 
with e = 0, Eex = 50A , E F = 100 A , and Uo = 0, where 
Tf refers to (1 + — B^). (c) The spin-down transmission 
coefficient T[ in a FIS junction with e = 0, E ex = 50Ao, 
E F = IOOAq, and U = 0, where 7\ refers to (1 + A L - B L ). 



sharply as we increase the barrier potential. For E ex < 
Ep, the increase of conductance with increasing x can 
not happen at x < 7r/2. Therefore, we would need a very 
large barrier potential to let x > 7r/2 in order to observe 
the first chacteristic of the Klein tunneling for E ex < 
Ep. This implies that obseving the first chacteristic of 
the Klein tunneling for E ex > Ep is easier. We also 
find that E ex can affect the amplitude of oscillation. For 
E ex < Ep, the amplitude of oscillation decreases as E ex 
increases. The amplitude of oscillation become zero as 
E ex = Ep. However, for E° x > E ex > Ep, the amplitude 
of oscillation increases as E ex increases. For E ex > E^ x , 
the amplitude of oscillation decreases with increasing E ex 
again. We could also alter the amplitude of oscillation 



increasing x, while for E ex < Ep, the conductance drops by applying Un [28l| . However, unlike E ex , Uq would only 
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decrease the amplitude of oscillation. As E c ex > E ex > 
Ep, the exchange energy would enhance the amplitude of 
oscillation. In addition, the phase shift of ir/2 is not seen 
by adjusting Uq. These interesting features of G(0)~x 
curve can be understood as follows. 



C. Transmission vs incident angle 

From Fig. 7(a), it can be seen that the transmission 
coefficient for a NIS junction is insensitive to \ as the 
incident angle becomes very small. Here, we regard the 
transport at very small angle as the nearly normal in- 
cidence. Also note that the transmission coefficient T 
refers to (1+A—B) for a superconducting junction, where 
A and B denote the Andreev and normal reflection coef- 
ficients, respectively. However, A = \a\ 2 cos(a')/ cos(a) 
and B = \b\ 2 for a NIS junction, where a' is the Andreev 
reflection angle and a is the incident angle. The trans- 
mission coefficient for a normal metal-insulator-normal 
metal junction (NIN) also has the same feature[3, @|. 
Fig. 7(b) and (c) show that the transmission coefficient 
also depends on the spin of incident electrons for the FIS 
junction. This is different from a NIS junction. In ad- 
dition, comparing Fig. 7(b) with Fig. 7(a), we find the 
critical angle changes under the effect of E ex . It is clear 
from Fig. 7(b) and Fig. 7(c) that the insensitivity of the 
transmission coefficient to x f° r nearly normal incidence 
exists also in the FIS junction, and this is not affected 
by the exchange spliting. We conclude that the barrier 
potential has little effect on the transport of Dirac-likc 
quasiparticles for the nearly normal incidence. 

In addition, Fig. 7(b) shows that the transmission 
sharply drops to zero as a is close to the critical an- 
gle regardless of x- The transmission is independent of 
X if the incident angle equates to or a c . The bar- 
rier strength x could affect the transmission only in the 
range of a = — a c . Furthermore, we find that T{x = 0) 
is larger than T(x = tt/2) at any incident angle for 
E ex < Ep. Conversely, for E ex > Ep, T(x) is always 
less than T(x = tt/2). For brevity, we do not represent 
the transmission of the case with E ex > Ep. This indi- 
cates that the contribution of any incident angle would 
increases the amplitude of oscillation (G max — G m i n ). In 
short, the increase of critical angle could enhance the 
amplitude of oscillation. 

However, as stated in Sec. II. A, we do not consider 
the evanescent wave functions when we calculate the con- 
ductance of the FIS junction. As stated above, in the low 



1.5 



bias regime, for E ex < E£ x , a'j 



From Eq. (8), we 



see that if E ex < Ep, decreases with increasing E ex 
and conversely, increases with increasing E ex . How- 



ever, for E% > E ex > E, a 



a 



\3- 



a^ 3 . If the critical angle is smaller, the conductance is 



less sensitive to x an d the amplitude of oscillation would 
decrease. Therefore, in Fig. 6, the G(0)~x curve would 
display such behavior (as stated above) under the effect 
of E P „- 
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FIG. 8: (Color online) The conductance spectra (a) for several 
values of E ex with Ef = 0.5Aq and Uo = 500Ao and (b) for 
several values of Ef with E ex = 0.5Ao and Uo = 500Ao. 
The curves for \ = and \ = ir are identical. Note that Go 
is not a constant and £ = tvHvf/Ao. As e < (E ex — Ef), 
Go = A(E ex /A )(e 2 W/h^ ) and as e > (E ex - E F ), G = 
4((e + £F)/A )(e 2 W7/i£o), where ^ = nhv F /A . 



Now, we discuss the case of E ex = Ep = 100A be- 
cause it is interesting that the spin polarization is 1.0. 
The critical angle aS equates to sin~ 1 (\e + E ex — Ep\/\e + 



E ex would reduce E 



f -f E ex \). For the small bias voltage (e -C Ep), 
a c = s\Tr l {\e + E ex -E F \/\e + E F + E ex \ w 0. The Dirac 
fermion only proceeds with nearly normal incidence and 
hence the transmission coefficient is almost independent 
of x- Moreover, the conductance for E ex = Ep is also 
almost independent of x- 
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D. Conductance in the small Fermi energy limit 

Now let us discuss the case with comparable Ep,E ex 
and Aq. Here, we should take a large Uq to fulfil 
the requirement of the mean field approximation. For 
Uo » E F ,A , we let sin(7) — ► in Eq. (A9) and Eq. 
(A10). We find that a a and b a are independent of \- The 
conductance also becomes independent of x- From Fig. 
8, we find that the conductance spectra have singular 
points at eV = \Ep — E ex \, E ex + Ep, and Ao. Among 
them, the singular point at eV = Ao can be found in 
the conductance spectra of all kinds of superconducting 
junctions. Therefore, we do not discuss it any more and 
focus on the other singular points instead. The singular 
points at eV = \Ep — E ex \, E ex + Ep are just the cross- 
ing points of the conduction and valence bands, i.e., the 
Dirac points shown in Fig. 2. The appearance of these 
singular points indicates that the transport properties of 
quasiparticles drastically change when they go through 
these Dirac points. In fact, at these Dirac points, the 
transitions between the specular Andreev reflection and 
Andreev retroreflection happen, as mentioned in Sec. II. 
C. The transitions result in the singular behavior of the 
conductance near these two singular points. The posi- 
tions of the singular points are related to the exhangc 
energy and Fermi energy. This indicates that by ob- 
serving the Andreev reflection spectrum and labeling the 
positions of the singular points, we can measure the ex- 
change energy and Fermi energy of graphenc. 



V. CONCLUSIONS 

In conclusion, we have investigated the transport prop- 
erties of a graphenc FIS junction within the BTK for- 
malism. Unlike the previous works on grephene FS 
junctions [HI, [36|, among other things, we study the 
spin-polarization of tunneling current at the TS inter- 
face and also invesigate how the exchange energy influ- 
ences the characteristic conductance oscillation behavior 
in graphene junctions. We find that the spin polariza- 
tion P of DOS in a graphene ferromagnetic segment is 
not always positively correlated with the exchange en- 
ergy and can be tuned all the way up to full polarization 
P = 1 (i.e., half-metallic state) with the applied voltage, 
and that graphene FIS junction can exhibit the unique 
negative tunneling spin polarization (Pp < 0). In ad- 
dition, we have calculated the tunneling conductance of 
the graphenc FIS junction in the thin-barrier limit. We 
shows in this limit that the tunneling conductance oscil- 
lates as a function of barrier strength \ with a period 
of 7r, and the exchange energy E ex can affect both the 
amplitude and phase of oscillation. For E ex < Ep, the 
amplitude of oscillation decreases as E ex increases. How- 
ever, for E% x > E ex > Ep, the amplitude of oscillation 
increases as E ex increases, where E£ x = 2Ep + Uo- For 
E ex > E^ x , the amplitude of oscillation decreases with 
E ex again. The curve for E ex > Ep is phase shifted 



by 7r/2 in x, compared with the curve for E ex < Ep. 
The conductance spectra always have singular points at 
eV = \Ep — E ex \, E ex + Ep and A . This property sug- 
gests that one could measure the exchange energy and 
Fermi energy of graphene by observing the Andreev re- 
flection spectrum and locating the positions of the singu- 
lar points. The unique oscillation originates from the rel- 
ativity of quasiparticles in graphenc. The wave functions 
of relativistic quasiparticles with a linear energy disper- 
sion relation do not decay in the high barrier. E ex can af- 
fect the amplitude of oscillation because the critical angle 
a£ varies with E ex . The transitions between retroreflec- 
tion and specular reflection happen at eV = \Ep — E ex \ 
and eV = E ex + Ep. The difference in transport prop- 
erties between specular Andreev reflection and Andreev 
retroreflection results in the singular behavior of the con- 
ductance near these singular points. 
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APPENDIX A: DERIVATION OF SPIN- 
POLARIZED-DIRAC-BOGOLUBOV-DE-GENNES 
EQUATION 

The graphene Hamiltonian including both valleys K 
and K' is given by 

H± = —ihvp(<T x d x ± (Jydy) + U. (Al) 

Taking exchange splitting into consideration and intro- 
ducing real spin degrees of freedom, the graphenc Hamil- 
tonian can be rewritten as 

H' = H®l 2x2 + hU xi ®S z . (A2) 

H' acts on an eight-dimensional wave function with two 
sublattices, two valleys and two real spin components, 

* = 5> a ®&, (A3) 

cr 

with 

ct=(;),ci=(;). (A4) 
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The motion of the quasiparticlcs in a superconductor 
should be described by the Bogolubov-de Gennes equa- 
tion 



H' — Epl SxS 
A* 



Al8x8 

-{T'H'T- 1 -E F i 8x8 ) 



•q> h 



(A5) 



Here, T" is the time-reversal operator. When the spin 
degrees of freedom are also considered, the time-reversal 
operator reads 



r = 



-iSy)C, 



(A6) 



S z and S y are Pauli matrices acting on real spin space 
and —iS y reverses real spin directions. It can be shown 
that T'H'T'- 1 = H <£> i 2X 2 - hUxi ® S z . The solu- 
tion of Eq. (A5) must satisfy the condition: * /l = T^ e 



Therefore, by writing "J 6 = (u^ 

(vt v 



T U T u l u i 



T U T "1 u i 



') and *' 1 = 

), we find that the Bogolubov-de Gennes 



equation can be decomposed into four decoupled sets 



H± — Epl 2 x2 — Pahl2x2 

A*l2x2 — {Hi 



A\ 2 x2 
Ef^2x2 + P<jh\2x2) 



PhS = 



(e - E F + p a E ex ) 

Hvf 



cos(a' s ). 



(A8) 



Here, q is the transverse wave vector of the particle and 
we have assumed that all transverse wave vectors are 
equal, e is the excitation energy and ± denote the trav- 
eling directions of the particles. In other words, and 
travel the +x direction, while ipea arm V^Ts- travel the 
—x direction. a a is the incident angle of the electron, 
while a' s is the Andreev reflection angle of the hole. 

In order to obtain the wave functions in the insulating 
region, we simply replace Ep with Ep — Vq and let E ex = 
in Eq. (Al) and Eq. (A2). Therefore, in region I, the 
wave functions are given by 



4>t = (l,±e d 



\0 7 0) T e ±tp ' x e l<lv , 



= (0,0,l,Te ±ie ') T e ±lphX e lqy , 



(A9) 



with 



sin(0) 



sin(0') 



hvpq 



(e + Ep- Vo) 



hvpq 



(e-E F + V ) ' 



(A7) 



. (e + E F - Vp) 

Pe = Z COS(0), 

hvp 



This is the Eq. (3) in Sec. II A. 



APPENDIX B: CALCULATION OF THE 
ANDREEV AND NORMAL REFLECTION 
AMPLITUDES 

By solving the DBdG equation with ferromagnetism 
[Eq. (3)], we find the wave functions in each region of 
the model FIS junction shown in Fig. 1. In region F, the 
wave functions are given by 

= (l,±e ±M -,0,0) T e ±lp " : V™, 



= (0,0, l,Te ±W *) T e ±lphsX e tqv , (A8) 



with 



sin(a (T ) 



hvpq 



(e + Ep + Pa E ex ) 



. (e - E F + Vp) a 

Ph = z cos(0 ). 

nvp 



(A10) 



The indices a and a disappear because our insulator is 
nonmagnetic. 6 (6') denotes the incident angle of the 
electron (hole) in region I. The meanings of the other 
symbols are the same as that in region F. 
In region S, the wave functions are given by 

_ ^ e ±*/3 _|_ e ±i(7+/3) e —i(j> ± e i(±~f~4>)^T e ±ip 3 x-Kx e iqy 

(All) 



with 



. cos 1 (s/A ) if e < A Q 
P 1 -icosh _1 (e/Ao) if e > A 



sin(7) 



hvpq 
Ep + W 



sin(a4) 



hvpq 



(e - E F + p a E ex ) ' 



(Ep + Uo) . . 

Ps = ^ C0S( 7 ), 



(e + E F + p a E ex ) 

Pea = z cos(a CT ), 

hvp 



(E F + U Q )A Q . 



(A12) 
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Both ipt an d tyj travel in the +x direction. -0+ repre- 
sents the state traveling in the same direction as its wave 
vector while if>~ means the state traveling in the direction 
opposite to its wave vector. 7 is the angle of incidence of 
the quasiparticlcs in region S. 

Then, the overall wave functions in region F, region I, 
and region S read 



(A13) 



Note that a„ and b a are the amplitudes of the Andreev 
reflection and normal reflection, respectively, as spin-cr 
electrons incident. They can be used to calculate the 
tunneling conductance and we can obtain them by ap- 
plying the boundary conditions, 



^Fa\x = -d = */L=-d, */|x=0 = * 



S\x=0- 



(A14) 



In the thin-barrier limit, we let d — > and Vq — > 00 and 
introduce the finite barrier strength \ — Vod/hvF such 
that 9 — » 0, & — > 0, k e d — > — x, and k e d — ► x- Then, a CT 
and b„ can be written as 



APPENDIX C: CALCULATION OF THE 
TUNNELING SPIN POLARIZATION 



In order to obtain the tunneling spin polarization, we 
must calculate the tunneling spin current for the inci- 
dent spin-cr electron first. For a s-wave superconductor, 
a cooper pair is composed of two electrons with opposite 
spin directions. The Andreev reflected hole has the same 
spin direction er as the incident electron but travels along 
the —x direction opposite to the incident electron, as de- 
cribed in Sec. II. A. The motion of Andreev reflected 
hole with spin-cr is equivalent to that an electron with a 
tunnels through the F-I-S junction interface. Therefore, 
the Andreev reflected hole carries the spin-cr current {I s ) 
while the incident and normal reflected electron carries 
the spin-cr current {!„). Therefore, the difference and 
sum of spin current for the incident spin-cr electron are 
given, respectively, by 



de / da^Go- (e) cos(cv)x 
^0 



(1 - M £ ,a CT )| 2 ^^ ~ M^«-)| 2 )> 
cos(cv) 



2cos(o; CT )cos(7) _ l2 _ l4> 
a a = e x e v , 



2cos(a cr )A 1(T - X a 

x a 



-i2 X 



(A15) 



Here, X a , X\ c are given by 

X„ = X\ a e 1 " + X 2 a, 

X la = [cos(/3) cos(7) + sin(/3) sin(7) cos(2x)] 
+i[sin(/3) - sin(/3) sin( 7 ) sin(2x)]e- M * , 
A 2cr = [cos(/3) cos(7) - sin(/?) sin(7) cos(2x)]e _M ' i 



(1 + K( £ ,a.)| 2 ^ 
cos(ao-) 



da a G a {e) cos(acr)x 
\b a {e,a a )\ 2 ). (A17) 



Then, summing up the contributions of the incident spin- 
up and spin-down electrons, we obtain the tunneling spin 
polarization 



Pt = 



(Ais) 



+z[sin(/3) + sin(/3) sin(7) sin(2x)]. 



(A16) 
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